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Abstract. We demonstrate that the general second-order scalar-tensor theories, which 
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acceleration, could be strongly constrained from the argument of the gravitational Cherenkov 
radiation. To this end, we consider the purely kinetic coupled gravity and the extended 
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1 Introduction 

Cosmological observations of type la supernovae [1-3], the cosmic microwave background 
[4, 5], and the large scale structures [6-9] indicate that the universe undergoes a phase of 
accelerated expansion, and this discovery opened up a new field in cosmology. A number of 
attempts to explain the origin of the present cosmic acceleration have been proposed over 
the past decade. The Einstein's cosmological constant might be a possible solution, however, 
the smallness of value of the cosmological constant can not be explained naturally [10, 11]. 

The possible solution for an accelerated expansion of the universe at the present time is 
an alternative theory of gravity. So far various modified gravity models have been proposed 
such as the scalar-tensor theories [12-16], f(R) gravity [17-20], Dvali-Gabadazde-Porrati 
(DGP) braneworld model [21, 22], and Galileon gravity [23-32]. In these models, additional 
degrees of freedom can mimic the cosmological constant and lead to cosmic acceleration today. 
Most of these theories are a subclass of the most general second-order scalar-tensor theory, 
which was first constructed by Horndeski [33] and also independently derived by Deffayet 
et al. [34] as an extension of galileon theory. The most general second-order scalar-tensor 
theory is applied to the late-time accelerated expansion [35] as well as the inflationary models 
[36-39]. 

The Lagrangian in the most general second-order scalar-tensor theory contains the cou- 
pling of the scalar field 4> and its kinetic term X = —g^ u 'SJ tl (j)SJ v (f)/2 with gravity, such as 
G^cj), X)R and G^{4>, X)G flu V fl V / (j), where G<i((f>, X) and G^{(f), X) are arbitrary functions of 
<j) and X. This theory is covariant, but in the presence of a cosmological background Lorentz 
invariance could be broken due to these coupling terms. As a result, the propagation speed 
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of gravitational waves differs from the speed of light and also depends on the cosmological 
background. 

When the propagation speed of gravitational waves is less than the speed of light, the 
gravitons could be emitted through a similar process to the Cherenkov radiation [40-42] . The 
observation of the high energy cosmic rays puts constraints on this process, i.e., the speed 
of the gravitational waves. Assuming a galactic origin for the high energy cosmic rays, the 
lower bound on the propagation speed of gravity from gravitational Cherenkov radiation is 
given by [41, 42] 



where ct is the propagation speed of gravity. When the origin of the high energy cosmic rays 
is located at a cosmological distance, the constraint is four orders of magnitude tighter than 



In the present paper, we show that the argument of the gravitational Cherenkov radi- 
ation puts a tight constraint on general second-order scalar-tensor models on a cosmological 
background with a time-varying propagation speed of gravitational waves. As a demonstra- 
tion, we consider two models: the purely kinetic coupled gravity [43] and the extended galileon 
model [44], which are a subclass of the most general second-order scalar-tensor theory. 

This paper is organized as follows. In section 2, we briefly review the most general 
second-order scalar-tensor theory and the tensor perturbations. In section 3, we derive the 
gravitational Cherenkov radiation on a cosmological background. In section 4, we explic- 
itly show that gravitational Cherenkov radiation reject the purely kinetic coupled gravity 
model. In section 5, we briefly review the extended galileon model and see how gravitational 
Cherenkov radiation can tightly constrain model parameters. Section 6 is devoted to con- 
clusions. In appendix A-C, we summarize the scalar perturbations, derived in [36], and the 
coefficients of the scalar and tensor perturbations in various regimes in the extended galileon 
model, derived in [44]. In appendix D, a useful constraint on parameters in the extended 
galileon model is demonstrated. 

Throughout the paper, we use units in which the speed of light and the Planck constant 
are unity, c = h = 1, and Mp\ is the reduced Planck mass related with Newton's constant by 
Mpi = l/v / 87rG. We follow the metric signature convention ( — ,+,+, +). 

2 The most general second-order scalar-tensor theory 

The most general second-order scalar-tensor theory is described by the action, 



c-c T <2x KT 1& c 



(1.1) 



(1.1). 
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where K, G%, G4, and G5 are arbitrary functions of the scalar field (f) and the kinetic term 
X = — g^V ^cfN u 4>/2, Giff, and Gix stands for dGi/dcp and dGi/dX, respectively, and C m is 
the matter Lagrangian. We assume that matter is minimally coupled to gravity. Note that 
for the case, G 4 = M|j/2, the La grangian £4 reproduces the Einstein-Hilbert term. 

We consider the tensor perturbations in the most general second-order scalar-tensor 
theory on a cosmological background, and briefly review the results in derived in [36]. We 
briefly review the tensor perturbations in the most general second-order scalar-tensor theory, 
derived in [36]. The quadratic action for the tensor perturbations can be written as 



where 



c(2) _ 

F T = 2 
Gt = 2 



1 



dtd 3 xa 3 



G T h% 



Ga — X yfiGsx + 
G4 — 2XG4X — X ( H(pG$x — 



Here an overdot denotes differentiation with respect to t, and H 
eter. We find the propagation speed of the tensor perturbations, 
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Ft 
Gt 



(2.3) 

(2.4) 
(2.5) 

a/ a is the hubble param- 
(2.6) 



When G4 = G^cj)) and G5 = 0, the propagation speed of gravitational waves is equal to the 
speed of light. On the other hand, the propagation speed of gravitational waves depends 
on the cosmological background in the presence of G5 or G4 being dependent on X. If 
the propagation speed of gravitational waves is less than the speed of light, it is tightly 
constrained from gravitational Cherenkov radiation. 



3 Gravitational Cherenkov radiation in an expanding universe 

In this section, we derive the gravitational Cherenkov radiation in a cosmological background. 
For simplicity, we consider a complex scalar field with the action 

S m = J d 4 x^=g~ [-g^d^d^ - m 2 ^*^ - . (3.1) 

Here we assume the conformal coupling with spacetime curvature £ = 1/6, for simplicity, but 
this term can be neglected as long as we focus on the subhorizon scales, p/a,m 3> H, where 
p is the comoving momentum. The free part of \E' can be quantized as 



#(77, x) 



1 



d 3 p 
a J (2vr) 3 / 2 



(3.2) 



where 77 is the conformal time, b p and c P are the annihilation and creation operators of the 
particle and anti-particle, respectively, which satisfy the commutation relations [£>p,SL] = 
8(p ~ p')> [cp, cL] = <5(p — p'), and the mode function obeys 



(J~2 +P 2 + m 2 a 2 ^j rp p (rj) = 0. 



(3.3) 
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The WKB approximate solution is given by (e.g., [45]) 



cxp 



% \ %{rl)dri 



'/iH 



(3.4) 



with Q p (r)) = \/p 2 + m 2 a 2 . The WKB approximation is valid for 



1 d 2 n n 3 1 fdCl 



fi p dr] 2 2nj\dri 



(3.5) 



which can be satisfied as long as p/a,m 3> H. 

On the other hand, the action of the graviton is given by eq. (2.3), then, we have the 
quantized graviton field 



h 



\1V 



a V ~Cfr V 



d 3 k 

(2vr) 3 / 2 



e$a k h k ( V )e^ + e$athUri)e 



-ikx 



(3.6) 



where e$ is the polarization tensor, and ak are the creation and annihilation operators, 
which satisfy the commutation relation [cik,a k /] = S(k — k'), and the mode function satisfies 

d 2 



dr]' 



+ c 2 k 2 --)h k ( V )=0. 



For the case c s ~ 0(1) and c s k/a 3> H, we may write 



exp 



i / u k (r]')dr]' 



(3.7) 
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where we defined uj^ = c s k, and the approximate solution is valid as long as c s k/a 3> H. The 
interaction part of the action (3.1) is given by 



Si 



dtd i xah i jd i ^dj^ 



I 



dr/d 3 x hij dixfjdjifj* , 



(3.9) 



(3.10) 



In order to evaluate the gravitational Cherenkov radiation, we adopt the method de- 
veloped in [46, 47]. Based on the in-in formalism [48], the lowest order contribution is given 
by 



where we defined tp = a^>, and the interaction Hamiltonian is 

Hi = a [ d 3 xh ij d i ^d j ^*. 



(Q(t)) = i 2 r dt 2 r dh ([#/(*!), [iT/(* 2 ),Q]]) . 

Jt[ n •> tin 



(3.11) 



We consider the expectation value of the number operator and the initial state with the one 
particle state with the initial momentum, i.e., 6p in |0). Then the lowest-order contribution of 
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Figure 1. Fcynman diagram for the process 



the process so that one graviton with the momentum k is emitted from the massive particle 
with the initial momentum pi n , as shown in fig. 1, is written as [49] 

4 (A) 4 A) ) = 2K / dt 2 [ 2 dh (ffi(*i)4 (A) 4 A) #/(*2)) • (3.i2) 

J ti n ''tin 

Then, the total radiation energy from the scalar particle can be estimated as E = Y2\ Sk( a; fc/ a ) 
(4 «k /' wn ^ cn l ea( is to 

d?k ujk 



(2vr) 3 a 



v I 2 

dm -7-T \l— h k(Vl )% } (Vl ) V>p in (m ) eij'PinP/ 



a{vi) V £t 



(3.13) 



where pj + k = p; n (p^ + /c* = p- n ). With the use of the relation ^aN^iiF/I = Pi n snl ®i 
we have 



E 



d 3 k 



(2vr) 3 a 



w fc 4 ■ 4 

p in sm i 



\ 1 

dm 



a(m) V £r 



(3.14) 



We are now interested in the subhorizon scales, k/a, p/a, m, c s k/a S> H, and the situation 
so that the scale factor a is constant, then we can approximate as 



dm tztt a / w- h k (mWpf (m (m ) 



Vin la {mnQr 

i nr i i 
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Then the total radiation energy eq. (3.14) reduces to 

_ 1 /• (i 3 A; pf n sin 4 g27rr 
~~ 4£ T a 3 7 (2vr) 3 fyf^ a~~ 

Here we assumed the long time duration of the integration, 



dr)i exp [i(Q iQ - fi/ - w*;)(r/i - r/ ini )] . (3.15) 



(5(Q in - Qf - w fc ), 



(3.16) 



dr]i exp [i(O in - 12/ - u k )(m ~ ^ni)] 
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2vrT 



5(0^-0/-^), (3.17) 



- 5 - 



where T/a = rj — rj m . Then, we have the expression in the relativistic limit of the massive 
particle, p\ n /a S> m, 



dE p 2 n f°° dkk 2 rl 



J —— ^cos^sin 4 ^^-^-^). (3.18) 
dt 4y T a Jo 27T J-i 

Now consider the delta-function, which can be written as 

S (n in - n f - Uk ) = 2n f 5(n 2 f - (o in - u k f) (3.19) 



where uik = c s k, Q- m = ypf n + a 2 m 2 , and fl/- = a/ (pin — k) 2 + a 2 m 2 . With the use of the 
fact 

n 2 - (n in -u k ) 2 = -2 Pin k (cos 0-j- (1 ~ c J )fc ) , (3.20) 



where we defined (3 = p m /y P m + m 2 a 2 and p? n = |pi n | 2 , we find (cf. eq.(3.2) in reference by 
Moore and Nelson [41]) 

dE pi f k -dkk s . n , e 



dt iG T a 4 Jo 2vr 
with 



"^1 + ^ < 3 ' 22 ' 



and 

= ±-c 2 V 

Assuming (3 ~ 1, we have fc max ~ 2p- m / (1 + c s ) and 



7^2 fl-l) 



dt SttQto- 4 
which yields (cf.[41]) 



2 /-fcmax / k \ 2 

4(1 - c s ) 2 / dkk ( 1 - -— , (3.24) 



rfg GWpf n 4(1 - c s ) 2 
dt a 4 3(1 + c s ) 2 ' 



(3.25) 



where we introduce the Newtonian gravity constant by Gn = 1/16ttGt- One may notice that 
this definition of the Newton's constant is slightly different from that in the most general 
second-order scalar-tensor theory (cf. [50]), however, it does not affect the constraints signif- 
icantly. Our results are consistent with those in Ref. [41]. Then, a particle with momentum 
p cannot possibly have been traveling for longer than 

a 4 (1+c,) 2 1 
1 ~ G~ N W^f- (3 ' 26) 
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Therefore, the highest energy cosmic ray put the constraint on the sound speed of the graviton 
(,_*)< ^Ur^"^)*^)* (3 . 27) 

Since we are considering the theory on a cosmological background, the sound speed 
of the graviton is determined by the cosmological evolution of the background field. This 
situation is slightly different from that in ref. [41]. However, as we have shown in this section, 
the theory on a cosmological background can be constrained from the gravitational Cherenkov 
radiation when the speed of the graviton is smaller than that of light. Also, there are no 
the higher order nonlinear interaction terms of the graviton like the galileon cubic term that 
becomes important at short distance [51], which suggests that the nonlinear interactions of 
the gravitons can be ignored. 



4 Purely kinetic coupled gravity 

We first consider the modified gravity model, whose action contains a nonminimal derivative 
coupling to gravity. The action proposed by Gubitosi and Linder [43] is given by 



S = <Fx 



R + X + —G^V^Vut, 
pi 



(4.1) 



where A is a dimensionless constant. In this model, the arbitrary functions in eq. (2.2) 
correspond to K = X, G3 = 0, G4 = M|j/2, and G 5 = — Xcp/Mpy Using the matter 
density parameter O m = pm/SM^H 2 , the modified Friedmann equation can be written as 



1 = Q m + fL., where 



% = T77#777 (1 + 18C) • (4.2) 



Here we defined the key parameter, 



2 



A H 

C =1T72-> C *> ( 4 - 3 ) 

iVipj 

where C* = —1/18. The second inequality is the condition which ensures the positivity of the 
energy density of the scalar field, $7^ > 0. Using the gravity equations and the energy density 
P0 and the pressure for the scalar field, the effective equation of state, w e g = p^/ p^, can 
be written as 

= 1 + 30C 

WcS ~ 1 + (24- 6fy)C + 108(1 + n 4> )C 2 ' ( ' j 

Gubitosi and Linder showed that if the deviation parameter at the present time, 5 = (C* — 
C)/C* | 2= o, satisfies 5 < 2/5, corresponding to the condition for negative pressure w c g < 0, 
the kinetic term X behaves as the cosmological constant around the present time. 
The propagation speed of gravitational waves (2.6) can be written as 



2 _ Mpj + 2\X/M 2 l 
CT ~ M 2 x -2\X/M 2 { 
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The condition for avoiding ghosts of the tensor perturbations, Qt > 0, is 5 > 0,^(^1^ — 3), 
which is automatically satisfied, while the condition for avoiding instability c|n > is 

<4 - 6) 

Therefore, 5 > is required for avoiding ghost-instability. Thus the theoretically allowed 
parameter range is 

< 5 < \, (4.7) 
5 



which is equivalent with 



1<C(; = 0)<-1. (4i 



The propagation speed of gravitational waves in terms of £7^ is rephrased as 

The constraints from gravitational Cherenkov radiation ct > 1 — e, where e = 2 x 10 -15 , 
reads 5 > 1 — 0(e) from eq. (4.9), which contradicts with the condition (4.7). Equivalently, 
from eqs. (4.3) and (4.8), A is always negative, therefore, the propagation speed of gravita- 
tional waves is always smaller than unity from eq. (4.5). Thus this purely kinetic coupled 
gravity is inconsistent with the constraint from the gravitational Cherenkov radiation for any 
theoretically allowed parameter A. 



5 Extended galileon model 

In this section, we consider the model proposed by De Felice and Tsujikawa [44], which is an 
extension of the covariant galileon model [52]. In this model, the arbitrary functions has the 
following form, 

K = -c 2 M 2 4(1 - p2) Xf 2 , 
G 3 = C3 M 3 1 - 4p3 X P3 , 

G 4 = \Ml-c±M 2 -^XV\ 

G 5 = 3c 5 M- (1+4p5) Xf 5 , (5.1) 

where q and pi are the model parameters and Mj are constants with dimensions of mass. 
We impose the conditions that the tracker solution is characterized by H<p 2q = const and 
the energy density of the scalar field is proportional to 4> 2p . These conditions enable us to 
reduce the model parameters, which is given by p2 = p, P3 = p + (2q — l)/2, p4 = p + 2q, 
and P5 = p + (Qq — l)/2 1 . Note that the covariant Galileon model corresponds to p = 1 and 
q = 1/2. 

1 Kimura and Yamamoto considered the case : p = 1, q = n — 1/2, C4 = 0, and C5 = [53]. 
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5.1 Cosmological Dynamics 

In this subsection, we briefly review the background dynamics in the extended galileon model. 
For convenience, we write the mass dimension constants as 



M 2 
M 3 

M 4 

M 5 



(ff dS M P1 ) 1 /2 ; 

Mpi l-2p 3 \ V(1-4P3) 



#ds 2pa 

Mpi 2-2p4\V(2-4 P 4) 



#ds 2p4 

tt 2+2p 5 \ V(l+4ps) 
-HdS \ 



Mpi 1_2p6 



(5.2) 



where H^s is the hubble parameter at the de-Sitter point. At the de Sitter point H = and 
<j> = 0, we obtain the following relations from the gravitational and scalar field equations 



C2 



C3 



3(3a -4/3 + 2) / 2 



' as 



V2[3(p + q)(a- f3)+p] 
2p + q-l 



MS 



p+<? 



(5.3) 



where x = 4>/HM p \ and 



a 



(3 = 2V2 



4(2p4-l) ^Js x; 

2 

P5+1/2 



C4, 



MS 



C5- 



(5.4) 



Thus this model is characterized by only four parameters p, q, a, and /3. In order to simplify 
the analysis, we introduce the following variables, 



n 



(Xds\ 2q ( HdS 



\ x J 



H 



l+2q 



T2 



.1 y 1 

^dS / r\ 



p+29 

l + 2q 



(5.5) 



and the radiation density parameter Q r = p r /3H 2 M^. Note that the de Sitter fixed point 
corresponds to (ri,r2,O r ) = (1,1,0). 

Along the tracker n = 1, the evolution of r 2 and J7 r are governed by the following 
differential equations, 



, (l + s)(O r + 3-3r 2 ) 
r 2 = — — ! r 2 , 



9' 



sr 2 + 1 
£l r — 1 — 3r 2 — 4sr 2 
sr 2 + 1 



(5.6) 
(5.7) 



where a prime denotes a derivative with respect to = In a and only one parameter s = pjlq 
determines the background dynamics in the case of the tracker solution. In this case, the 
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density parameter of the scalar field is simply given by £1$ = r 2 , satisfying the constraint 
1 = Q<j) + VL m + Q, r . Integrating these equations yields the following algebraic equations, 

r 2 = b ia * {1+s) n r 1+s , (5.8) 
b ia 4{1+s) n r 1+s = 1 - fi r (l - b 2 a) , (5.9) 

where the integration constants are given by 

1 — n m Q — £l r o £l m o 
h = ^j+ s , b 2 = -—, (5-10) 

and £l m o and Q r o are the matter and radiation density parameter at present, respectively. 
To see how the Friedmann equation is modified, we rewrite the algebraic equation (5.9) in 
terms of the hubble parameter H, then we find 

H \ 2 / H \ ~ 2s 

—J = (i-n m0 -n r0 ) ( —J +n m0 a- 3 + n r0 a- 4 . (5.11) 

This modified Friedmann equation is known as the Dvali- Turner model [54]. The authors in 
[53] placed the observational constraints on this modified Friedmann equation (5.11) in the 
special case p = 1 using type la supernovae and the CMB shift parameter and showed that 
the model parameter s has to be small, s <C 1, in order to be consistent with cosmological 
observations 2 . 

5.2 Conditions 

In this subsection, we summarize the theoretically allowed parameter space in the extended 
galileon model, discussed in [44] , and show that the constraint from gravitational Cherenkov 
radiation is crucial. To avoid ghost-instabilities, we must impose the conditions, Qt > 0, 
4 > 0, Qs > 0, and c| > in the history of the universe. The coefficients in the tensor 
and scalar perturbation equations in terms of r%, r2 r , and the model parameters are listed 
in appendix B. We find that the propagation speed of gravitational waves along the tracker 
r\ = 1 is written 



2 _ 2(1 -2p- Aq){2q+pr 2 ) + 3a(2g + pr 2 )r 2 - 3/3(1 - 2p - 4g)(3 - 3r 2 + q.)r 2 
° T ~ (1 - 2p - 4q)[2 + 3(a - 2(3)r 2 ](2q + pr 2 ) 



(5.12) 



Note that eq. (5.12) reduces c\ = 1 when a = (3 = 0, which correspond to G4 = Mpj2 
and G5 = 0. We further impose no-instability condition at r 2 = r2 im i n , where a minimum 
of propagation speed of gravitational waves is located. Setting T\ = 1 and Q r ~ 0, the 
minimum of c 2 ^ is given by eq.(5.12) at r 2 = r 2 ^ m i n , 



r 2. 



2(3 + 2p)(l -2p- 4q)q (3 - 8pq(p + 2q)a ± V3T7 



/r 2 , (5.13) 



2 Observational constraints on eq. (5.11) from type la supernovae, cosmic microwave background, and 
baryon acoustic oscillations including the cosmic curvature K in the context of the extended galileon model 
has been recently studied by De Felice and Tsujikawa [55]. They found that the parameter s is constrained 
to be s = 0.0341°;^ (95% CL) in the flat case K = 0. 
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where 



ri = (1 - 2p - 4q) (p + 2q)q (3 x [4(p + 2q) (p -3qa)a 

+2(1 -2p- 4g){(3 + 2p) - 3(3 - 4g)/3}/3 + 3[3 - 16g(l - 2q) - 2p(3 - 8q)]a 0\, 
T 2 = 4p 2 (p + 2q)a - 18(p + 2g)(l - 2p - 4g)/5 2 + (1 - 2p - 4g)[2p(3 + 2p) + 9(p + 2q)a]p. 

(5.14) 

The conditions for avoiding ghost-instabilities in the regimes along the tracker are given by 

£7s|ri=l,r 2 «l > 0, Gs\de Sitter > 0, 

c l'|ri=l,r2<l — 0' c llde Sitter > 0, 

^T|ri=l,r 2 <l > 0, Gride Sitter > 0, 

2 i 2 i 

Op\ ri =l,r2<.l — 0) c Tlde Sitter > 0, 

c T|r2,min > 0. (5.15) 

If the initial condition of T\ is r\ <C 1, we then must impose the conditions for avoiding 
ghost-instabilities in the regime r\ <C 1 and f2<l, which is given by 

£s|ri-Cl,r 2 -Cl > 0, c|| ri< i jr2< i > 0, 

^r|ri<l,r 2 <l > 0, CT\rx<t.l,T2<.l — 0- (5.16) 

We also impose the condition that the other fixed points r a and r& (see appendix C) is not 
real or outside the interval < r\ < 1, which is given by 

A < or r aib < or r a<b > 1. (5.17) 

Note that as long as the initial condition of r\ is near r\ = 1 and the scalar field follows the 
tracker from early stage, these conditions (5.16) and (5.17) do not have to be imposed. 

Let us classify the constraints into four classes: (a) the constraint from the gravitational 
Cherenkov radiation, which is given by eq. (1.1) and eq. (5.12) with setting ct = ct\ z =o, (b) 
the theoretical constraint (5.15) to avoid the ghost-instabilities when the scalar field follows 
the tracker solution from early stage, assuming that the tracker is near r± = 1 initially, 
(c) the theoretical constraint (5.16) and (5.17) in addition to (5.15) to avoid the ghost- 
instabilities when the scalar field does not follow the tracker solution initially, assuming that 
the initial condition of n is sufficiently small, (d) the other constraint from the cosmological 
observations, type la supernovae, the shift parameter from the cosmic microwave background, 
and the baryon acoustic oscillations. 

Figure 2 shows the allowed regions to satisfy the constraint (a) and the constraint 
(c) for p = 1 and q = 1/2 (left panel) and p = 1 and 5/2 (right panel), where we adopt 
n mQ h 2 = 0.1344 and n r0 h 2 = 4.17 x 1CT 5 with h = 0.7. In this case, we see that there is 
no overlap region except for a = and f3 = 0. Thus, the constraint from the gravitational 
Cherenkov radiation is crucial. Figure 3 is the same as figure 2, but for the constraint (a) and 
the constraint (b). We see that the allowed region in parameter space is significantly reduced, 
by combining with the constraint from gravitational Cherenkov radiation (a). Especially, 
there is no overlap region with the positive values of a and j3, in figure 3. In general one can 
show that both the constraints (a) and (b) impose a and (3 to be negative or zero for any 
values of p > 1 and q > (see appendix D). 
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Figure 2. The allowed parameter space which satisfies the constraint (a) from the gravitational 
Chcrcnkov radiation (1.1) and the constraint (c) from (5.15), (5.16), and (5.17). The left panel 
assumes p = 1 and q = 1/2, while the right panel does p = 1 and q = 5/2. 




u 1 z -4 -2 2 4 6 

a a 



Figure 3. The allowed parameter space which satisfies the constraint (a) from the gravitational 
Chercnkov radiation (1.1) and the constraint (b) from (5.15). The left panel assumes p = 1 and 
q = 1/2, while the right panel does p = 1 and q = 5/2. 

We must further include the constraint from cosmological observations (d) . The authors 
in ref. [56] investigated the constraint on the covariant galileon model (p = 1 and q = 1/2) 
from the observational data of type la supernovae, the shift parameter from the cosmic mi- 
crowave background, and the baryon acoustic oscillations. They showed the early tracking 
solution, corresponding to the case of (b), is disfavored by the cosmological constraint (d). 
On the other hand the solutions that approach the tracker solution only at late times, cor- 
responding to the case of (c), are favored (also see [55]) taking small spatial curvature into 
account. However, the latter case is significantly constrained by combining the constraint 
(a), though we do not take the spatial curvature into account. 

Thus, the constraint from the gravitational Cherenkov radiation plays a very important 
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role to reduce the allowed parameter-space of the extended galileon model. In ref. [57], it 
is demonstrated that the integrated Sachs Wolfe effect derives a stringent constraint on a 
subclass of the galileon model. Further tight constraint could be obtained by combining these 
constraints. 

6 Conclusion 

In this paper, we studied constraints on the general scalar-tensor theories on a cosmological 
background, whose propagation speed of gravitational waves differs from the speed of light, 
using the survival of high energy cosmic ray against the gravitational Cherenkov radiation. 
In these theories, the coupling of the scalar field eft and its kinetic term X with gravity 
causes the violation of Lorentz invariance in a cosmological background, leading to a time- 
dependent propagation speed of gravitational waves. We demonstrated that such a model 
can be constrained using the survival of high energy cosmic ray against the gravitational 
Cherenkov radiation. 

We first considered constraints on the purely kinetic coupled gravity and found that 
the conditions for the existence of a desired late-time solution and avoiding ghost-instability 
is0<<5<2/5 while the constraint from the gravitational Cherenkov radiation gives 5 > 
1 — 0(e), where e = 2 x 10 -15 . Thus the purely kinetic coupled gravity is inconsistent with 
the argument of the gravitational Cherenkov radiation. 

We also focused our investigation on the extended galileon model, which is a generaliza- 
tion of the covariant galileon model in the framework of the most general second-order scalar- 
tensor theory. We showed that there is no allowed parameter space except for a = f3 = by 
combining the condition for avoiding ghost-instabilities and the constraints from the grav- 
itational Cherenkov radiation if the initial condition of r\ is sufficiently small. Even if the 
initial condition of r\ is placed near the tracker r\ = 1, the allowed parameter space is tightly 
constrained by combining the gravitational Cherenkov radiation and cosmological constraint 
such as type la supernovae, the shift parameter from cosmic microwave background, baryon 
acoustic oscillations. 

Thus the constraint from the gravitational Cherenkov radiation is important to con- 
strain the general second-order scalar-tensor theories on a cosmological background, whose 
propagation speed of gravitational waves is less than the speed of light. 
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A Scalar perturbations 

Here we summarize the scalar perturbations including coefficients in the most general second- 
order scalar-tensor theory derived in [36]. For the unitary gauge 4> = </>(*) with the line 
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element ds 2 = -N 2 dt 2 + "f i j(dx i + iW£)(cfa^' + Wdt), where N = 1 + a, JVj = <9;/3, and 
7ij = a 2 (t)e 2 ^5ij, the quadratic action for the scalar perturbations after solving the constraint 
equations can be written as 



4 2) = / dtd 3 xa 3 



QsC 2 - §(VC) 2 



a- 



(A.l) 



where 



Ss = J^St + 3Ct, (A.3) 

and 

6 = -<j)XG 3X + 2HG 4 - 8HXG 4X - 8HX 2 G 4XX + <pG 4<t> + 2X<j)G^ x 

-H 2 <p {5XG 5X + 2X 2 G 5XX ) + 2HX (3G 5(f> + 2XG 5(f)X ) , (A.4) 

S = XK X + 2X 2 K XX + l2Hj>XG 3X + 6H<pX 2 G 3XX - 2XG 3<t> - 2X 2 G 3<t>x 
-6H 2 G 4 + 6[H 2 (7XG AX + 16X 2 G 4XX + 4X 3 G 4XXX ) 
-Hj>(G H + 5XG 4(j>x + 2X 2 G^ XX )]+30H 3 ^XG 5X + 26H 3 <j)X 2 G 5XX 

+4H 3 <j ) X 3 G 5XXX - 6H 2 X(6G 5<t> + 9XG 5( p X + 2X 2 G 5(j)XX ). (A.5) 

The propagation speed of the scalar perturbations is defined as 

4 - (A.6) 

B Coefficients and propagation speed in various regimes 

In this appendix, we summarize the coefficients and propagation speed in the tensor and 
scalar perturbation equations in the extended galileon model in various regimes, derived in 
[44]. 

In the regime, r\ = 1 and r 2 <C 1, the coefficients (2.5), (2.6), (A.3), and (A.6) are given 

by 

Gs\ ri =l,r 2 <l ^6q[p~ 3(q - 2fi)q] r 2 , (B.l) 
c|| ri=1>r2<<1 ~ {4 P 3 (Q r + 3) - 2p 2 {(n r + 3)(6/3 - 3a + 2) 

-2q[3Q r + 11 - 3(a - 2/3)(ft r + 3)]} - 3{/3(fi r + 3) 

+8q 3 (n r + 5)(a - 2/3) - 2< ? 2 (7ft r + 27) (a - 2/3) 

+g[3a(fi r + 3) - 2/3(5fi r + 17)]} 

-p{(n r + 3) (3a - 12/3 - 1) + 4g 2 [(a - 2/3)(9ft r + 33) - 2(ft r + 5)] 
+9(12(20 - a)(30 r + 10) + 6fi r + 22]}} 

xl/[24g 2 (2p + Aq - l){p - 3(a - 2/3)g}], (B.2) 

g T |r 1= i,r 2 «i ~ ^ [2 + 3(a - 2/3)r 2 ] , (B.3) 

4k=i,r- 2 <i - 1 " {6[2(a - 2/3)^ + 3/3]p + 24(a - 2/3)c/ 2 

+ 3/3(16? - 3) + 3/3(2p + 4g - l)O r }/4g(2p + 4q - l)r 2 . (B.4) 
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At the de Sitter point, n = r 2 = 1, the coefficients (2.5), (2.6), (A. 3), and (A. 6) are given by 

_ 6(p + 2g) (3a - 6/3 + 2) [p - 3(a - 2/3)g] 
^ |deSittcr " [2p - 6(a - 2/3)g - 3a + 6/3 - 2]2 ' l±4 - Dj 

4lde Sitter = {6/3 + 4p 2 + p [9(a - 2/3) 2 + 3a - 12/3 + 4g(6/3 - 3a + 2) - 2] 
+3(a - 2/3)[3/3 + q(9a - 12/3 - 8q + 6)]} 
3(2/3-a)(2< ? + l) + 2p-2 



6(6/3 - 3a - 2)(p + 2g)(2p + 4? - l)(p - 3a<? + 6/3g) ' 



(B.6) 



StU Sitter = i (3a -6/3 + 2), (B.7) 
2 | 2(2p + 4g-l)-3a 

C T |de Sitter - ( ^ + 4? _ ^ _ 6/3 + 2 ) " ^ 

In the regime, n <C 1 and r 2 <C 1, the coefficients (2.5), (2.6), (A. 3), and (A. 6) are given by 

Ssk«i,r 2 «i * 3(p + 3<z)(2p + 6g - l)/3rf- 1)/(2g+1) r 2 , (B.9) 
p + 3(7 — 2 

" 2(p + 3 g )(2p + 6 g -l) (1 + ^ ' (B - 10) 

aT|n«l,r 2 <! * 1 - 3/3r 2 ri p - 1)/(29+1) , (B.ll) 

2 I _ i , 3(4p + 12g - 5 - 3» r ) (p-i)/(2 9 +i) rR 19 x 

crlri<i,r 2 «i — 1 H 4 p + i2 g - 2 1 (B.12) 

C Other fixed points 

There also exist the other fixed points found by De Felice and Tsujikawa [44], which is 
characterized by the equation, 

p(3a - 4/3 + 2)r 2 + [2/3(p + 3q) - 3a(p + 2q)]n + 2/3(p + 3q) = 0, (C.l) 

where = r a and r&, and 

_ 3a(p + 2q)-2/3(p + 3g)± 
ra ' b ~ 2p(3a -4/3 + 2) ' 1 > 

A = [2/3(p + 3?) - 3a(p + 2q)} 2 - 8/3p(3a - 4/3 + 2)(p + 3g). (C.3) 

D Constraint on the values of a and /3 

In this appendix, we assume p > 1, q > 0, and r 2 = 1 — f2 m o — f^rCb where fl m o/i 2 = 0.1344 
and ri r o/i 2 = 4.17 x 10 -5 with h = 0.7. One can show that a and /3 must be negative or zero 
to satisfy both the constraints (a) and (b) for any values of p > 1 and q > 0. This can be 
proved as follows. The upper bound of the constraint (b) is determined by the straight line 
in the plane of a and /3, 

B = -a - f - 1 r , (D.l) 
H 2 3(2? + 1)' V ; 

which comes from c||d e sitter > (see e.g., the left panel of figure 3). On the other hand, the 
constraint (a) is characterized by the two straight lines, 
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and 



a = 2{p + 2q)(2q+pr 2 ) 

P (2p + 4q-l){4q-3 + 3r 2 + 2pr 2 -n r ) ' 1 ' ' 
Note that the lines of eqs. (D.2) and (D.l) are parallel each other. Since we have 
2{p+2q)(2q+pr 2 ) 1 

(2p + 4q-l)(4:q-3 + 3r 2 + 2pr 2 -n r ) 2' 1 '' 

at present epoch, the lines of eqs. (D.2) and (D.3) intersect at a point with a > and (3 > 0. 
(D.4) means the slope of the line (D.3) is larger than that of the line (D.l), then the lines of 
eqs. (D.l) and (D.3) intersect at a point with a < and (3 < 0, Therefore, a and /3 must be 
negative or zero to satisfy the constraints (a) and (b) in the case of p > 1 and q > 0. 
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